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We investigate the explicit renormalization group for fermionic field theoretic representation of 
two-dimensional random bond Ising model with long-range correlated disorder. We show that a 
new fixed point appears by introducing a long-range correlated disorder. Such as the one has been 
observed in previous works for the bosonic (y)*) description. We have calculated the correlation 
length exponent and the anomalous scaling dimension of fermionic fields at this fixed point. Our 
results are in agreement with the extended Harris criterion derived by Weinrib and Halperin. 



Critical properties of systems with short-range and long-range correlated randomness have been studied extensively 
[1-8]. One important question to address is whether the introduction of weak randomness changes the universality 
class of transition. That is, in the renormalization group (RG) language, whether the disorder is relevant at the critical 
point of pure system or not. According to the well known Harris criterion [3,6], disorder is irrelevant if dv > 2, where 
d is the dimensionality and v is the correlation length exponent of pure system. This criterion should be modified 
in the presence of long-range correlations in the disorder. A special type of such a disorder has been considered by 
Weinrib and Halperin. They showed that the disorder with power law correlation ~ x^^^'^ ( for large separations x ) 
is irrelevant if [{d — 2p)v — 2] > for p > 0, whereas the usual Harris criterion recovers for p < [7]. Therefore the 
existence of long-range correlations in the disorder would have significant effect in the sense that they can change the 
universality class of phase transitions. 

The most useful method to study disordered systems in RG language is the Replica method. By employing this 
method, one can average out the free energy using a trick based on the equation \nZ = lim„_^o ^ ' '^^^ idea is 
then to average Z". However, RG analysis can be implemented explicitly without averaging on disorder [14]. Using 
this method, one can avoid some of the mathematical problems e.g. the n — > limit in the replica method. Moreover, 
application of the explicit method in some cases, such as the one studied in this letter, is more straightforward. 

In this paper we consider random-bond Ising model with fermionic action and long-range correlated disorder. The 
effect of short-range correlated disorder has been studied in previous works [9-14]. Also the ip^ version of this model 
is investigated through the double expansion near four dimensions with short-range and long-rang disorder [7,8]. 
Wc found a new long-range fixed point and we calculate both correlation length exponent and scaling dimension of 
fermionic field at this fixed point. 



Two dimensional Ising model near its critical point can be described by the massive free fermionic action 



I. INTRODUCTION 



II. EXPLICIT RENORMALIZATION GROUP 




(1) 



X 



is a two component Grassmannian field (-01 = ■00 and 



(2) 
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where are Pauli matrices. The two point function of the Grassmannian fields of the model in the momentum space 



IS 



Go{p){27TrS{p + q) =< mm >= -^^{27r)H{p + q), (3) 



where '>p{p) is the Fourier transform of ip{x) 

ip{p) = ilj{x)exp{—ip ■ x) (4) 



and 



fp_ 



'^(^) = / / 77^i^ip)exp{ip-x). (5) 



The randomness can be inserted into the action in the following way 

X 

where c(x) is a random variable with the following correlation function in the momentum space 

< c{p) >= 

< c{p)c{q) >= [Do + Dplpr^nSip + q). (7) 

Here Do and Dp are short-range and long-range disorder strengths respectively. Positivity of two point function 
imposes some restrictions on Dq and Dp for example Dp should be positive for p > and Dq should be positive for 
p < 0. Here we consider the case with < p < 1. 



^=\( ( -^ZT?^(-P)(^P + ^)^iP) + t; f f f f 4^4^^(Pi)c(-pi -P2)^/'(P2), (8) 



We can write the action in the momentum space as 
where 

c{p) = c{x)exp{—ip.x). 

X 

It is convenient to introduce a diagrammatic representation. We have the following vertex 



Pi +P2 



Pi P2 



Wavy line stands for ipj solid line for ip and zigzag represents the c insertion. Momentum conservation should be 
regarded at each vertex. So we have the following graphs for vertex renormalization 



^AAAA/^ ^AAAA/^ ^ 4^24: ^A/W^/^ - 



where the dashed line represents tjjip propagator. Symmetry factors will be canceled with nl of perturbation expansion 
in each order. The key equation will be [14] 



c'{ri,r2) ^ c{ri + r2) - c{ri - q)Go{q)c{q + r2) + / / c(ri - 9i)Go(gi)c(gi - g2)Go('?2)c((72 + ^2) 

J q J qi •> q2 

c{ri - qi)Go{qi)c{qi - ^2)60 (g2)c(g2 - 93)60(93)2(93 + r2)- (9) 



qi ■'q2 -J qs 

The random function in the original action has zero mean and we should keep the mean value fixed after RG 
transformation. So the mean value of new function should be extracted and in fact it will renormalize the mass and 
the kinetic term. Then we need to introduce a field renormalization constant Z to keep the coefficient of kinetic 
term equal to i, just as in the original action. After rescaling momentums (r —>■ Ar) and fields ("tp — > ^^/^V")? 
renormalized random function and mass will be 

Cr{puP2) = -^[C'{pi/\P2/X)~ < C'{p,/X,p2/X) >]. (10) 

rUr ^ ^{m+ {c')r=o) (11) 

A simple expansion of the second term of (9) in powers of r, with to = 0, leads to the following expression for field 
renormalization constant to first order in disorder strength 

Z = ^ + ^i>^''-^) + 0{D'). (12) 

The short-range correlated disorder does not contribute up to this order. 
So the renormalized mass is 



TTlr 



Ato |l - ^(In A)i?o - ^^i>^'' 1)^P + ((^)' - {h + h+ h)DoD, - [h + h)Dl + •••}, (13) 



where the /j's are 



h 
h 

h 
h 
h 



qi Jq: 



-92 -(92 - 9i) + m'^ 



qi •Jq2 



{q( +rn'^){q^ + 'm?){{q2 - qiY + nn?) 

{-qi.{q2-qi)+m^)\qi\-^P 
(qf + m'^){q'^ + m'^){{q2 - qiY + m^) 

(-91-92 +™^)|92 -9ir^'' 



qi Jq2 ill + "^^)(9i + m'^){{q2 - 9i)^ + m"^) 

(-9i.(92+gi)+m^)|gi|-^P|g2-9ir'^ 
(qj + m?){ql + ■m?){{q2 - qi)^ + m?) 
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(-9i-(9i + 2^2) +TO^)|92 - 91 1 91!"^" 



qi J q2 







> 




91 


- 92I 


A' 






> 


TT 


9i 


-92I 


A' 






> 


TT 


91 


- 92I 


A' 






> 


TT 


l9i 


- 92I 


A' 



(gf +?7i2)2(q| + m2) 



we did not include the contribution of Z in the second order. 

The renormalized values of disorder strengths can also be obtained by calculating the correlations of renormalized 
random function Cr{pi,P2)- Up to second order of bare strengthes, the renormalized strengthes are 



InA 1 + p 

In A 1 + p 



Dor = D,{1 —Do 2^^^" 1)^" + ■••)' (14) 



D,.. = X'^D,{1 —Do ^^i>^"' l)Dp + •••). (15) 

By differentiating these equations with respect to In A, and then replacing the bare parameters in terms of renormalized 
ones, one can obtain the following Wilson's functions 



^=™,(l-i^I,.„-i^if,.+---). (16) 

Si + + .-.). (17) 

dD„r ,11 + , , , 

^,-D,.i2p^-D..^-^D,. + ...). (18) 

It is clear from the equations that at the Gaussian fixed point {Dq = D* = 0), Dq is marginaUy irrelevant and 
Dp is relevant. So, by introduction of small amount of long-range correlated disorder it becomes unstable in Dp 
direction. Apart from the trivial Gaussian fixed point, we see that there is a nontrivial fixed point at Dq = and 
Dp = jT^TT^"'"^''. The new fixed point is attractive in all directions in the Dq and Dp plane. The RG fiows starting in 
the vicinity of Gaussian fixed point, end up at the nontrivial fixed point. One of the most important critical exponents 
is the correlation length exponent, which turns out to be different at two fixed points, i/ = 1 at the Gaussian fixed 
point and v = at the nontrivial fixed point. The result is in agreement with [7]. A good numerical confirmation 
of this relation for the special choice of p = ^ can be found in [16]. 

Also the field renormalization would acquire an anomalous dimension which is defined through the asymptotic 
behavior of vertex function in the long wavelength limit as 

r(2)(p) ^pi-'j (19) 

so at the nontrivial fixed point we have 

d2__^^ (20) 

Now we want to show that, some less singular terms in the correlation functions of c^, which were omitted in the 
equation (7), are irrelevant. We start with a general form of the correlation function of disorder 

< c{p) >= 0, 

n 

<c{p)c{q) >= (Do + Y,DpM'''')S{p + q), (21) 
1=1 

and then find the Wilson's functions for the Dq and Dp. 

1=1 

'^-DpA2p.-'-n,.r-±i±^Dp^,. + ...). (23) 

i=i 

Here we have n nontrivial fixed points at Dq — and all D*. — except one of them, say D*^ which is j^^tt^^^''' . 
At these fixed points the RG eigenvalues are —2pi and 2{pi — pi) which means that there is just one stable fixed point 
(with pi > Pi for all i). Other fixed points are unstable at least in one direction. 

Finally we want to compare the above results with the replica results. In the case of short-range disorder, as pointed 
out by Murthy [14], the results of the replica method and the explicit method are in complete agreement. The replica 
action for the model considered here can be easily obtained, 

5 = ^ ^ 1p{x){d + m)'llj{x) - ^ X! ^a{x)i'a{x)i^l3{x)iJl3{x) - 
x,a x,OL,f3 

x,y,a,l3 ^' 

The third term in this action is nonlocal, and the computation of the (3 functions are not easily tractable. This 
shows the advantage of explicit calculations with which one can avoid such a nonlocal action. The other advantage 
of this method is that there is no restriction on the distribution function of disorder while it should have a Gaussian 
distribution in the replica method. 
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